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Abstract 

We derive constraint on the effective number of neutrino species Ny from the cosmic microwave 
background power spectrum of the WMAP and galaxy clustering power spectrum of the SDSS 
luminous red galaxies (LRGs). Using these two latest data sets of CMB and galaxy clustering 
alone, we obtain the limit 0.9 < Ny < 8.2 (95% C.L.) for the power-law ACDM flat universe, with 
no external prior. The lower limit corresponds to the lower bound on the reheating temperature 
of the universe Tr > 2 MeV. 
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I. INTRODUCTION 



The standard model of cosmology with the concordance set of parameters can suc- 
cessfully reproduce a broad range of the cosmological data such as the big bang nucle- 
osynthesis (BBN), the cosmic microwave background (CMB) anisotropics and large scale 
structure (LSS). The relativistic degrees of freedom present after the BBN epoch in the 
standard cosmology are photons and three generations of neutrinos. In models of particle 
physics/cosmology, however, there are many candidates that could additionally contribute 
to the relativistic components of the universe: sterile neutrinos IL gravitational waves, 
(pseudo-)Nambu-Goldstone bosons such as axions and majorons ^, and the active neu- 
trinos themselves if they have large lepton asymmetries ^. Furthermore, the energy density 
of the relativistic particles can be smaller than in the standard cosmology, if the thermal- 
ization of the neutrinos are ineffective as in the MeV-scale reheating scenarios In 
fact, in a certain class of models, especially those accompanied by the late-time entropy 
production j?], the (final) reheating temperature tends to be quite low, and it often falls in 
the MeV range. Therefore it is of great importance to study the possible effects of varying 
the effective number of the relativistic particles on the cosmological observations, not only 
to make the standard cosmology more established, but also to probe and constrain a certain 
class of models in the particle physics/cosmology. 

Recent precise observations of the CMB anisotropics and LSS make it possible to mea- 
sure the relativistic degree of freedom in the universe through its effects on the growth of 
cosmological perturbations. These effects come from the fact that the density perturba- 
tion does not grow (the gravitational potential decays) during the radiation-dominated era. 
Specifically, more relativistic degree of freedom causes more early integrated Sachs- Wolfe 
effect on the CMB power spectrum, which leads to higher first peak height. Also, since it 
delays the epoch of the matter-radiation equality and makes the horizon at that time larger, 
the turnover position of the matter power spectrum is shifted to larger scales and the power 
at smaller scales are suppressed. Therefore, by observing CMB and LSS, we can measure 
the relativistic degree of freedom during the structure formation. In detail, assuming the 
smallest scale relevant to our observations to be about 5 Mpc, since the structure formation 
of that scale begins around the temperature T ^ 20 eV (at which the scale enters the hori- 
zon), these observations probe the relativistic degree of freedom at T < 20 eV. Thus, CMB 
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and LSS can measure the relativistic degree of freedom independently of another well-known 
probe, BBN, which measures it in much earlier universe around T = O(MeV). 

In this paper, we analyze the most recent data sets of CMB and LSS, respectively using 



the Wilkinson Microwave Anisotropy Probe (WMAP) 3-year data j^, ^, 10, 11| and the 
Sloan Digital Sky Survey (SDSS) luminous red galaxies (LRGs) power spectrum data |12| . 
and would like to discuss the constraints from them. The WMAP data now can give clear 
features of the first and second peaks of the CMB power spectrum. In particular, the 
precision around the first peak has been already cosmic variance dominated and so has been 
the measurement of the early integrated Sachs- Wolfe effect. However, this is not the case for 
the measurement of the relativistic degree of freedom since it is almost completely degenerate 
with the value of the matter density. We would like to demonstrate this degeneracy by the 
WMAP data alone analysis. Then, as is well known, since the matter power spectrum 
has somewhat different degeneracy pattern from the CMB, it is broken by combining the 
CMB with the matter power spectrum data. We would like to see how and how much 



the degeneracy is broken by adding the matter power of the SDSS 



analyses of earlier data sets are found for example in Refs. [13 



data, in Refs. 
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19|, 



3, El using WMAP 1st year data, and 



RG sample. Related 
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23|, 



2J] using 



WMAP 3-year data, with various combinations of astrophysical data sets and priors for 
cosmological parameters. 

The new point of our analysis on the relativistic degree of freedom in the universe is to use 
the power spectrum of the SDSS LRGs 13]. Although there are similar analyses using the 
power spectrum of the SDSS main galaxies [3] and/or the one of the 2dF Galaxy Redshift 
Survey (2dFGRS) [26|, it is important to revisit the issue with the new power spectrum data. 
This is because not only have the LRGs more statistical constraining power (the effective 
volume of the LRG survey is about 6 times larger than that of the SDSS main galaxy sample 



and over 10 times larger than that of the 2dFGRS [l2:]), but also there seems to have 



Deen a 



tension between the power spectra of the 2dFGRS and the SDSS main galaxies [26|, 127]]. The 



discrepant measurements of the relativistic degree of freedom between these two samples as 



found in Refs. 



(we show the result of Ref. [22| in Table [TTl Compare the 4th and 



5th lines) are considered to be caused by this tension in the power spectra. It is shown 



in Ref. 



27l | that the discrepancy is due to the scale dependent bias which was not taken 



into account in the SDSS main galaxy analysis. The LRG analysis in Ref. 12] models this 
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effect of scale dependent bias in the same way as the 2dFGRS analysis and they found the 
extracted cosmological parameters, especially the matter density, are in excellent agreement 
with those from WMAP alone and WMAP + 2dFGRS. Therefore, it is useful and of great 
importance to investigate the relativistic degree of freedom using the LRG power spectrum 
and see whether the discrepancy in its derived value from the two galaxy surveys is resolved 
to give a reliable constraint.^ 

We describe our analysis method in Sec. [TTl and our results are presented in Sec. IIIII In 
Sec. llVt we briefly review a cosmological scenario with low (MeV-scale) reheating tempera- 
ture and explain how the relativistic degree of freedom is modified compared to the standard 
case. In Sec. El we discuss our result in relation to previous works and see how much the 
current observations allow the relativistic component of the universe to deviate from the 
value in the standard cosmology. We also emphasize its implication for the lower bound on 
the reheating temperature. 

II. ANALYSIS 

The quantity we try to constrain in this paper is the effective number of neutrino species, 
Ny, which is widely used to quantify the energy density of the relativistic component in 
the early universe. It is given by = (prei — P-y) / Pu,thm, where is the photon energy 
density, prei is the total energy density of photons, three active species of neutrinos and 
extra relativistic contribution, and pi,,thm is defined as p„^thm = (77r^/120)(4/ll)'^/^T^ using 
the photon temperature T after the electron-positron annihilation. Pu,thm corresponds to 
the energy density of a single species of neutrino assuming that neutrinos are completely 
decoupled from the electromagnetic plasma before the electron-positron annihilation takes 
place and they obey Fermi-Dirac distribution. 

We constrain in the fiat ACDM universe with the initial perturbation power spectrum 
which is adiabatic and described by power law. This model has 6 cosmological parameters, 
the baryon density Ub, the matter density the normalized Hubble constant h, the reion- 
ization optical depth r, the scalar spectral index of primordial perturbation power spectrum 

^ At present, the scale dependent bias is modeled in a very plienomcnological manner and more detailed 
modeling is considered to be required. However, this is the on-going issue in the community and beyond 
the scope of our paper. 
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Hs and its amplitude A {uo = Qh"^, where Q is the energy density normahzed by the critical 
density). Theoretical CMB and matter power spectra are calculate d by the CMBFAST code 
28| and by the likelihood codes of the WMAP 3-year data 0, H 11] and of the SDSS 



LRG powe. ~ data We apply ™deUng of no„-U„ea..ty and scale dependent 

bias as in Ref. [12] to the linear matter power spectrum before fitting to the LRG data. Since 
we omit the process of "dewiggling" (which is found to be justified in Sec. IIIip . this modeling 
has two parameters, galaxy bias factor b and non-linear correction factor Qni. Specifically, 
we connect the linear matter power spectrum P\in{k) and the galaxy power spectrum Pga\{k) 
by 

We calculate the functions of N,^ by marginalizing over the above parameters (6 

parameters for WMAP alone and 8 for WMAP+SDSS). The marginalization is carried out 
by the Brent minimization [2^ modified to be applicable to multi-dimension parameter space 
as described in Ref. [3 



III. RESULT 

We show the results of minimization in Fig. [H We give the values of some of the 
best fit cosmological parameters as functions of N^, in Fig. [2l We have checked that the 
results for standard three neutrino species agree with the WMAP js] and SDSS [3] groups' 
analyses. For the WMAP 3-year alone case, it has been checked in Ref. 311 that the best fit 
x"^ and parameters agree. With regard to WMAP and LRG combined analysis, our best fit 
parameter values for three neutrino species are Ub = 0.0222 ±0.0007, Um = 0.1288 ±0.0044, 
h = 0.718 ±0.018, r = 0.088 ±0.029, = 0.958 ±0.016, ag = 0.770 ±0.033 (we report here 
(Tg instead of A to compare with Ref. [13]), b = 1.877 ± 0.065 and Qni = 30.4 ± 3.5. The 
central values fall well within the la ranges of the constraints derived in Ref. 13] and the 
la errors are almost identical to those quoted in Ref. 1^. This empirically shows that the 
"dewiggling" we mentioned in Sec. [Ill for the non-linear modeling can safely be neglected as 
is documented in the likelihood code of Ref. [12] . This makes sense as follows. Since the 
process of the dewiggling mainly decreases the amplitude of the acoustic oscillations in the 
matter power spectrum, it mostly affects the parameter estimation of 00^. However, Ub is 
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FIG. 1: as functions of Ny. The red solid line uses the WMAP 3-year data alone and the 

green dashed line uses WMAP 3-year and SDSS LRG power spectrum. 

more precisely determined by the CMB when we do the combined analysis so neglecting the 
dewiggling does not affect much the parameter estimation using the present galaxy clustering 
data. 

The limits corresponding to A^^ = 4 are Ny < 25 for WMAP 3-year alone and 
0.8 < Ny < 8.0 for WMAP and SDSS LRG combined. Since functions show some 
asymmetric features, we derive 95% confidence limits by integrating the likelihood func- 
tions £ = exp(-AxV2)- This yields 95% C.L. bound of < 42 for WMAP alone and 
0.9 <Ny< 8.2 for WMAP+LRG. 

We observe that the CMB alone constraint is very weak and the LSS data significantly 
reduces the allowed region. We note that our WMAP 3-year limit is somewhat weaker than 
the earlier constraints quoted as CMB alone limit 

nni 

IigI . [isl . even though they 



use data before the WMAP 3-year release. We compiled them in Table [B For example, 



Ref. 



has derived A^^j, < 9 (95% C.L.) using the WMAP 1-year data alone, much more 



stringent than our bound < 42. We can ascribe this apparent discrepancy to the prior on 
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h adopted by Ref. [I6j, h < 0.9, which affects N^, constraint through the well-known h — 
degeneracy. We find this degeneracy in our analysis too as in Fig. [2] (b) (although we show 
the result for < 10, for example, N^, = 22 can fit the WMAP 3-year data with h ~ 1.3). 
They are positively correlated which is explained as follows. The effect of increasing on 
CMB power spectrum is cancelled by increasing Um so that the epoch of matter-radiation 
equality occurs at the same redshift. Then, h has to be increased so that Qm ~ 0.25 leading 
to acoustic peaks at observed positions in the fiat universe. These features are explicitly 
demonstrated by the solid lines in Fig. [2] (a)-(c). 

The degeneracy is broken by adding LSS information as is clearly seen in Fig. [H This can 
be understood by another well-known fact that the shape of the matter power spectrum is 
determined by the combination Qmh rather than uj^n = ^mh"^- When Um is varied, obviously 
it is impossible to find h which preserves both Qm and Qmh. Thus, when we try to fit the 
CMB and LSS simultaneously, h faces the dilemma of fitting the CMB (Qm) or the LSS 
{flmh). We can see this in the panels (c) and (d) of Fig. [2l 

Our final result, the constraint on Ni, from the WMAP 3-year CMB power spectrum and 
the SDSS LRG power spectrum, can be summarized as 

0.9 <K< 8.2 or = 3.1^11 (2) 

at 95% C.L., whose center value is quite close to the standard model of three active neutrino 
species. We will be discussing the constraint in connection with other works in Sec. |Vl 



TABLE L Summary of CMB alone limits. 





CMB data 


95% limit 


Prior on h 


Hannestad [14] 


pre- WMAP 


< 19 or 24 


OA<h< 0.9 


Bowen et al. [15] 


pre- WMAP 


0.04 <N^ < 13.37 


0.4 < /i< 0.95, h = 0.65 ± 0.2 Gaussian 


Crotty et al. [16] 


WMAPl 


N^<9 


0.5 <h< 0.9 


Hannestad [18] 


WMAPl 


< 8.8 


0.5 <h< 0.85 


Barger et al. [19] 


WMAPl 


0.9 <N^< 8.3 


0.64 < /i< 0.8 


This paper 


WMAP3 


< 42 


NONE 
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WMAP 
WMAP+LRG 




FIG. 2: The best fit values of some of the cosmological parameters as functions of A'^^. The red 
solid lines are for the WMAP 3-year data alone and the green dashed lines are for WMAP 3-year 
and SDSS LRG power spectrum combined. 

IV. EFFECTIVE NUMBER OF NEUTRINOS IN THE UNIVERSE WITH LOW 
REHEATING TEMPERATURE 

Before we move on to discuss our result in the next section, it will be useful to review the 
low (MeV-scale) reheating scenario. In this scenario, the effective number of neutrinos iVj, 
can deviate from the standard value. We would like to briefly explain this scenario and how 

nn 

Ny and the reheating temperature Tr are related. For more details, we refer to Refs. [5|, 16[ ■ 
The standard big bang model assumes that the universe was once dominated by thermal 
radiation composed of photons, electrons, neutrinos, and their antiparticles. The reheating 
temperature is the temperature at which the universe becomes such radiation dominated 
state and it is usually assumed to be so high that every particle species is in thermal 
equilibrium. In particular, neutrinos are considered to obey Fermi distribution. 
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FIG. 3: Taken from the calculation in Ref. (a) The relation between the effective neutrino 
number Ni, and the reheating temperature Tr. (b) The solid line shows the ^He abundance 
as a function of the reheating temperature Tr. The dashed line is calculated with Fermi 
distributed neutrinos with Ni, of the panel (a) (namely, only the change in the expansion rate due 
to the incomplete thermalization is taken into account). The baryon-to-photon ratio is fixed at 
= 5 X 10-1°. 



What if the reheating temperature is lower, say, several MeV? In contrast to electrons 
that are always (at least until the temperature drops below a few eV) in thermal contact 
with photons via electromagnetic forces, neutrinos interact with electrons and themselves 
only through the weak interaction. The decoupling temperature of the neutrinos should 
be around 3 MeV for the electron neutrinos and 5 MeV for the muon and tau neutrinos, 
respectively (the difference comes from the fact that the electron neutrinos have additional 
charged current interaction with electrons). Therefore the neutrinos might not be fully 
thermalized and lead to N^, < 3 if the reheating temperature is in the MeV range. 

In fact, the reheating temperature as low as a few MeV can be found in many cosmological 
scenarios. To avoid the overproduction of the unwanted relics such as the gravitinos, one 
needs to require the reheating temperature low enough ^. In extreme cases it may be in 
the MeV range. Further, the thermal history of the universe may not be so simple that 
the universe might have underwent several stages of the reheating, and the final reheating 



^ This is the case if the gravitinos are thermally produced 3^, |33| . On the other hand, when the gravitinos 
are non-thermally produced by inflaton deca^, lower reheating temperature leads to more gravitinos, 



making the gravitino-overproduction severer 



35 



36 



37| 
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temperature may be very low. For instance, late-time entropy production [3] is one of the 
plausible ways to solve problems associated with the unwanted relics, and the reheating 
temperature often falls in the MeV scale. 



In Ref. jo], we have calculated how much neutrinos are thermalized when Tr = O(MeV) 
and have derived the relation between Tr and which is shown in Fig. [3] (a). Specifically, 
we have solved numerically the momentum dependent Boltzmann equations for neutrino 
density matrix, fully taking account of neutrino oscillations. For later convenience, we also 
show the ^He abundance Yp in the MeV reheating scenario in Fig. [3] (b). It should be noted 
that Yp increases while Ni, decreases in this scenario. This is in contrast to the conventional 
non-standard Ny scenario where decreasing Ny accompanies decreasing Yp. The difference 
occurs as follows. Since the latter assumes the thermal (Fermi) distribution for neutrinos as 
in the standard cosmology, only the expansion rate is modified and particularly it has the 
neutron-proton conversion rate identical to the standard one. Meanwhile, since the MeV 
reheating scenario makes the neutrino distribution less thermalized one, the neutron-proton 
conversion rate is significantly modified in addition to the expansion rate. To elucidate the 
effect of the modified neutron-proton conversion rate, we draw the dashed line in Fig. [3] (b) 
which expresses (fictitious) Yp when we include only the change in the expansion rate. 

We can now convert our constraint on A'^i,, Eq. ([2]), into the lower bound on Tr using 
Fig. [3] (a): 

Tr>2 MeV. (3) 

We will discuss this CMB-I-LSS constraint on Tr, paying particular attention to the com- 
parison with BBN bound, in the next section. 



V. DISCUSSION 



We have shown that new data of the SDSS LRG power spectrum can considerably shrink 
the allowed region of N^, from the one obtained using the WMAP 3-year data alone by a 
factor of six. In terms of the extra relativistic particle species other than three species of 
active neutrinos, the LRG data reduces the upper limit by a factor ^ 7.5, from 39 to 5.2. 
Moreover, combining with the LRG data gives a finite lower limit on the effective neutrino 
number, A^j, > 0.9. This translates into the lower bound on the reheating temperature of 
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the universe, Tr> 2 MeV as described in Sec. HVl 

A comparison to Ref. 2^ who has reported the constraints using earher data sets is in 
order. They have provided constraints on using various combinations of cosmological 
data sets, which we compiled in Table [TTl We can summarize their finding that their Lyg 



data [38] and/or the galaxy clustering power spectrum data from the SDSS main sample 25| 



prefer Ni, > 3 at more than 95% confidence level whereas the 2dF galaxy power spectrum 



261] does not show such a non-standard feature. Our new constraint is quite similar to the 
latter, WMAP3+2dF(+supernovae) constraint of A^^ = 3.2+^^ (95% C.L.). This resuh is 
reasonable since the SDSS main galaxy power favors significantly higher value of Qm than 

;he 2dF power [8] but the SDSS LRG power gives Qm which is close to the 2dF value 



12l |. The robustness of the estimation of from the SDSS LRG clustering is thorough 



y 



tested by means of the power spectrum shape |39| and the baryon acoustic oscillations |40|. 
Since galaxy clustering basically measures the matter-radiation equality, this robustness is 
considered to be transferred to our estimation of Ni,. We can conclude that although the 
constraints from both galaxy surveys has converged with central values around the standard 
value of three, allowed regions are large enough to cover the constraints obtained with Lya 
forest data whose central values are around 5. We have to wait for more study on the 
Lya forest analysis and future C MB/ LSS observations (the PLANCK sensitivity for N^, is 



4l|) to see whether the present Lya data would hint for 



forecasted to be 0.2, see e.g. Ref. 
non-standard physics.^ 

A cosmological constraint on N^, can also be obtained from the primordial '^He abundance 
Yp. While ^He has logarithmic dependence on the baryon-to-phton ratio rj, the only parame- 
ter in the standard BBN, it is very sensitive to N^y since it modifies the expansion rate during 
the BBN period and shifts the epoch of the neutron-to-proton ratio freeze-out. The deu- 
terium, D, constrains these parameters almost in the opposite manner. It is quite sensitive 



to 1] but has only mild dependence on N,^. More details are found in e.g. Refs. 13|, ll9|, |43 |. 
Actually, the BBN bound is more conventional than the structure formation constraint but it 
has somewhat checkered history since it is very difficult to estimate systematic errors for de- 
riving the primordial abundance from ^He observations. Although the D abundance is often 



^ It may suggest that the effective number of neutrinos increases after BBN. In Ref. 
such scenario is feasible by decaying particles. 



421, it is shown that 
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TABLE II: Comparison of N^, constraints using various data set combinations. "All" refers to 
WMAP3 + other CMB + Lya + galaxy power spectrum (SDSS main sample + 2dF) + SDSS 



baryon acoustic oscillation (BAO) + Supernovae la (SN). See Ref. 



221 for details. 



Seljak et al. [22] 



This paper 



95% limit 



A^. 



4.^ 



-1.6 
-1.4 



6.0^:^ 



+2.1 
-1.7 

h2.3 
-3.2 



= 3.9 



= 5.2 



+2.1 



A. 



"J-J--2.2 



Data set 



All 
All + HST 
All - BAO 
All - Lya 
WMAP3+SN+SDSS(main) 
WMAP3+SN+2dF 
All-2dF-SDSS(main) 



WMAP3+SDSS(LRG) 



or example, recent studies have shown the importance of underlying stellar 
, . This leads to significant increase in Yp and enlarged errors. Therefore, 



considered to more robustly probe the primordial abundance, since it does not have much 
sensitivity on Ny as mentioned above, systematic errors for Ny estimation are dominated by 
those of ^He. 
absorption [4^, 

the most recent Ny constraints Ny = (68% CL.) f4S|] has a higher central value and 

larger errors than the earlier results. Nevertheless, the current BBN bound is significantly 
tighter than our WMAP+LRG bound and is completely covered by our bound (this is not 
the case for a so-called MeV reheating scenario. The significance of BBN and CMB/LSS 
is reversed in this scenario. We comment on it below). At this stage, we can say that 
the present CMB plus galaxy clustering data provides a complementary constraint to the 
BBN. Our analysis shows consistency between the constraints derived from totally different 
physical processes and at distant epochs providing a strong support for standard cosmology, 
but relatively large error bars still leave some room for non-standard physics. 

Lastly, let us comment on the implication of our results for MeV-scale reheating scenar- 
ios. Since the reheating temperature is an important but not yet known parameter that 
characterizes the early evolution of the universe, it is valuable to derive an observational 
constraint. As shown in Sec. IIVI our estimation of Ny, in particular the lower bound of 
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A*",^, provides us with the lower bound on the reheating temperature as Tr> 2 MeV, once 



we use the relation between N^, and the reheating temperature given in Ref. 



a]. We cau- 



tion that on the contrary to CMB+LSS bound, the lower bound on obtained from Yp 



such as the one of Ref. 



43| we quoted above cannot be taken at face value. This is es- 



sentially because < 3 in low-reheating scenario implies not only less radiation density 
but also less neutron-to-proton conversion rate, which greatly affect the "^He yields by BBN. 
The latter effect was not taken into account when deriving the bound on N^j in Ref. 



That is why it cannot 



Yp = 0.249 ± 0.009 



43 



3e applied to the MeV-scale reheating scenarios. It turns out that 



44l |, which does not reject relatively large value of Yp, does not give 



meaningful lower bound on Tr (see Fig. E] (b) ). To derive a lower bound on Tr from BBN 
data alone, one needs a concrete upper bound on the "^He abundance, which is difficult to 
obtain due to a possibly large systematic error. More detailed discussion is given in Ref. . 
At present, CMB-I-LSS do better job in setting lower bound on Tr. It is quite intriguing 
that we have obtained the lower bound Tr > 2 MeV, which is just before BBN begins, even 
without resort to the BBN data. This is because the decoupling of the weak interactions 
accidentally occurs immediately before the BBN epoch. Due to this coincidence, we were 
able to derive the concrete and tight bound on Tr. So far, it has been considered that the 
observations on the light-element abundances are indispensable to probe the BBN epoch. 
However, our results unambiguously show that one can extract informations on the universe 
at T = O(MeV) by CMB+LSS data, and that one doesn't have to rely on the observed 
light-element abundances, which may have large systematic errors. When the PLANCK 
data becomes available, the lower bound can be improved up to ~ 5MeV, the decoupling 
temperature of the muon and tau neutrinos. 
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